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ON THE ABSORPTION OF X-RAYS IN THE INTERSTELLAR MEDIUM 
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ABSTRACT 

We present an improved model for the absorption of X-rays in the interstellar medium (ISM) intended 
for use with data from future X-ray missions with larger effective areas and increased energy resolution 
such as Chandra and the X-ray Multiple Mirror mission, in the energy range > 100 eV. Compared to 
previous work, our formalism includes recent updates to the photoionization cross section and revised 
abundances of the interstellar medium, as well as a treatment of interstellar grains and the H2 molecule. 
We review the theoretical and observational motivations behind these updates and provide a subroutine 
for the X-ray spectral analysis program XSPEC that incorporates our model. 

Subject headings: Atomic data — ISM: abundances — dust, extinction — X-rays: ISM 



and explain our choices of elemental and molecular abun- 
dances (§2.3) as well as dust grain parameters (§2.4). We 
summarize the model, compare it to previous models, and 
point out its remaining shortcomings in Section 3. In the 
appendices, we derive the influence of dust on (Tism (Ap- 
pendix A) and describe the implementation of our model 
into XSPEC (Appendix B). 

The intended readership of this paper consists of users 
of the new missions, many of whom do not intend to be- 
come specialists on the ISM. We have therefore attempted 
to include sufficient information to make this paper self- 
contained, and have chosen a review-like style. 

2. MODELING THE X-RAY ABSORPTION IN THE ISM 

2.1. X-ray Absorptivity in the ISM 

The total photoionization cross section of the ISM, uism, 
is obtained by summing over the contributions of the as- 
trophysically relevant elements. Taking into account the 
phases of the ISM, tiisM can be written as 

CiSM — Cgas + Cmoloculos + f grains- (1) 

As is common for X-ray applications, we normalize ctism 
to the total hydrogen number density (i.e., in molecular, 
neutral, or ionized form), A^h, so that the observed X-ray 
spectrum of a source. Jobs, is given by 

lobsiE) = e--^«"(^) /source (i?), (2) 

where A'h is measured in atoms cm~^ and /source (£') is the 
X-ray spectrum as emitted by the X-ray source. We drop 
the explicit energy dependence of (Tism for the rest of the 
paper. 

The contribution of the gaseous, non-molecular cold and 
warm phases of the ISM to ctism is obtained by summing 

^When referring to the cross sections of BM92 we use their subroutines as implemented as the phabs model in XSPEC, version 10. Oz 
(Arnaud 1996). This implementation contains improvements compared to the original publication. Starting with XSPEC, version 11.0, other 
cross-sections are used in the computation of this model (see §2.2 below). 
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1. INTRODUCTION 

The precise knowledge of the modification of X-ray ra- 
diation in the interstellar medium (ISM) is of crucial im- 
portance for the understanding of X-ray spectra from cos- 
mic sources since the observed X-ray spectra have to be 
corrected for this modification before interpreting the ob- 
served data. The current state of the art in the computa- 
tion of the energy-dependent photoionization cross section 
of the ISM, (Tism, has been summarized by Balucihska- 
Church & McCammon (1992, BM92). These authors 
greatly improved^ the precision in the computation of ctism 
compared to previous implementations (Strom & Strom 
1961; Brown & Gould 1970; Fireman 1974; Ride & Walker 
1977; Morrison & McCammon 1983, MM83) by allowing 
to adjust the assumed abundances of the ISM and by us- 
ing polynomial fits to the photoionization cross sections of 
Henke et al. (1982). 

Although the BM92 model provides sufficient precision 
for the data analysis with current instruments, the advent 
of X-ray missions with instruments of high energy reso- 
lution and large effective areas, such as Chandra or the 
Newton X-ray Multiple Mirror mission (XMM-Newton), 
will result in large improvements in the precision of X-ray 
astronomical measurements. It is therefore necessary to 
have tools at hand that reflect the current understanding 
of the physics of the interaction of X-rays with the ISM at 
the percent level. 

In this paper we improve on the models of BM92 by 
incorporating recent improvements in the photoionization 
cross sections and by taking into account the physics of the 
interstellar dust and the molecular phase of the interstellar 
medium. Section 2 describes the formalism employed to 
compute Ctism- We describe the atomic physics data (§2.2) 



2 



Wilms, Allen, & McCray 



the photoionization cross sections of the individual atoms 
and ions in these phases, weighting their contributions by 
the abundances. Hence, 

o-gas = ^Az- az,i • (1 - I3z,i) ■ (Thi{Z, i) (3) 

where we use the notation pioneered by Ride & Walker 
(1977), in which Az = iV(Z)/iV(H) is the abundance in 
number of element Z with respect to hydrogen, az,i = 
N{Z,i)/N{Z) is the fraction of ions of element Z that are 
in ionization stage i, and (Jhi{Z, i) is the total photoioniza- 
tion cross section of element Z in ionization stage i. The 
depletion of the elements into grains is taken into account 
through the "depletion factor" 1 — j3z,i- 

For the molecular phase of the ISM, only molecular hy- 
drogen needs to be taken into account due to its large 
abundance (~ 20% of the hydrogen might be molecular, 
see §2.3). Therefore, 

CTmolcculos = ^H2fbf(H2). (4) 

For material in the grain phase of the ISM, several com- 
plications arise. Since the X-ray optical depth of a large 

grain can be much greater than one (sec Fig. 3), most 
of the X-ray absorption will occur on its surface. Conse- 
quently, the contribution of material "within" the grain 
to cTisM is small. Thus, the net effect of this shielding is 
that the presence of grains reduces the effective absorptiv- 
ity of the ISM relative to the absorptivity of a completely 
gaseous medium. As we show in Appendix A, the total op- 
tical depth for grains with a size distribution dngr(a)/da 
is given by 

'^grains — *^grains-^H 

= A^h4 f ^^^cTgeom (1 - exp (-(a)(iV))) da 

(5) 

In the derivation of equation (5) we assume that the grains 
are of chemically homogeneous composition with an av- 
erage photoabsorption cross section {a), and can be ap- 
proximated as spheres with radius a and average column 
density {N) (measured in atoms cm~^). Since the grains 
are partly transparent to X-rays, their absorption cross 
section is smaller than the grain geometrical cross section, 
Cgeom = 7ra^. Finally, the number of grains per hydrogen 
atom along the line of sight is ^g. We discuss the validity of 
these assumptions in Section 2.4 and give explicit formulae 
for all parameters of equation (5) in Appendix A. 

In the following sections we describe the data used for 
evaluating Eqs. (1) and (5) in greater detail. 

2.2. Atomic Physics 

Despite decades of effort, the knowledge of elemental 
photoionization cross sections in the X-ray regime is still 
uncertain. In many cases, therefore, one has to rely on 
either theoretical work or interpolation along the isoelec- 
tronic sequence to obtain usable cross sections, and then 
use compilations of experimental data to gauge the pre- 
cision of these theoretical computations. In what follows 
we will describe the basis of the photoionization cross sec- 
tions used here, focusing on the region E > 100 eV. The 
emphasis is on the choice of cross sections that are good for 
instruments with a resolution comparable to the CCDs on 



Chandra or XMM-Newton. Therefore, we will not make 
an attempt to include the resonances that are important 
close to the absorption edges that might well be important 
for the grating instruments aboard these satellites. 

For the computation of the photoionization cross section 
of H, we use the fitting formula of Band et al. (1990), which 
provides a very good approximation of the exact, analytic 
cross section but is much faster to evaluate. The max- 
imum deviation between the Band et al. (1990) formula 
and the analytical cross section (e.g., Bethe & Salpeter 
1957, eq. 71.7, divided by two) is ~ 1.5% at energies be- 
low 2keV. 

The cross section for helium is taken from the recent 
theoretical evaluation of Yan, Sadeghpour & Dalgarno 
(1998). Autoionization resonances occur at energies less 
than 0.08 keV, and the four strongest resonances are in- 
cluded here (as in BM92), using values from Oza (1986) 
and Fernlcy (1987). For energies between 0.1 keV and 
1 kcV, the helium cross section is smaller by a factor of 
~ 0.9 to ^ 0.6 than that assmncd in the subroutine used 
by BM92, which is based on theoretical work of Chen, 
Cooper & Brion (1991) and on experimental data by Marr 
& West (1976). The Yan, Sadeghpour & Dalgarno (1998) 
cross section, however, agrees with other more recent the- 
oretical computations (e.g., those of Verner & Yakovlev 
1995) and is in agreement with more recent experimental 
data (Yan, Sadeghpour & Dalgarno 1998, and references 
therein). It also has the correct E^^'^ asymptote required 
from non-relativistic quantum theory (Bethe & Salpeter 
1957). For implications of this smaller He cross section on 
studies of the ISM, see, e.g., Arabadjis & Bregman (1999). 

For all other elements we use the photoionization cross 
sections as originally published by Verner et al. (1993) and 
later updated by Verner & Yakovlev (1995), and which are 
used in XSPEC version 11.0. These authors fitted the re- 
sults from Hartrec-Dirac-Slater computations of Cbf for in- 
dividual subshells to a semi-empirical fitting formula. The 
major advantage of this semi-empirical fitting formula over 
previous fitting formulae (such as, e.g., polynomial fits) is 
that it requires only a small number of fitting parameters 
to express (Tbt to high accuracy while preserving the cor- 
rect asymptotical properties. A comparison between the 
experimental data and these cross sections is presented by 
Verner & Yakovlev (1995), who focus on the overall agree- 
ment throughout the electromagnetic spectrum. In addi- 
tion, we compared the total photoionization cross sections 
computed from the fitting formula with the compilation of 
data presented by Henke, Gullikson & Davis (1993) for the 
range from 0.1 keV to lOkeV. Except for energies close to 
the K- and L-shell edges, our comparisons indicate max- 
imum relative deviations between these cross sections on 
the order of 5% above 1 keV. 

In their determination of Ubf close to the K- and L-shell 
edges, Henke, Gullikson & Davis (1993) used a linear inter- 
polation through the resonances to the edge energy. Since 
this edge energy was obtained from experimental data, 
it includes the effects of energy shifts in the solid state 
and thus is not representative of the true edge energy for 
the gaseous state of the materials. Indeed, the difference 
in edge energy between the solid and the gaseous state 
can be up to several 10 eV (see, e.g., Nicolosi, Jannitti 
& Tondello 1991, for an illustrative experimental exam- 
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pie). In many cases, Henkc, GuUikson & Davis also do not 
consider the different subshcU threshold energies. On the 
other hand, the Verner & Yakovlev (1995) cross sections 
are based on quantum-mechanical computations that in- 
clude the resonances, and the edge energy is unaffected by 
solid state effects. We have checked the precision of the 
Verner & Yakovlev edge energies by comparing them with 
the computations of Gould & Jung (1991) which are more 
appropriate for our work than the Henke et al. (1982) or 
Henke, GuUikson & Davis (1993) values since they assume 
the elements to be in the gas phase. For Z > 8 the edge 
energies of Verner & Yakovlev (1995) agree on the 10 eV 
level with those of Gould & Jung (1991). For C, N, and 
O, the edge energies adopted by Verner & Yakovlev (1995) 
arc between those of Henkc, GuUikson & Davis (1993) and 
Gould & Jung (1991). We therefore decided to adopt the 
Verner & Yakovlev edge energies for this work. Note, how- 
ever, that these differences in the edge energy arc relevant 
only for instruments with extremely high resolution and 
are irrelevant for most other work. 

For molecular hydrogen we adopt the cross sections re- 
ported by Yan, Sadcghpour & Dalgarno (1998) for energies 
above 85 cV. These cross sections arc claimed to be accu- 
rate to about 5% and have the correct E~^-^ asymptotics. 
Albeit outside of our formal energy range {E > 100 eV) 
we note that the fit formula given by Yan, Sadeghpour 
& Dalgarno (1998, eq. 18) for energies below 85 eV has 
significant deviations with respect to the tabulated cross 
sections of Samson & Haddad (1994), on which the Yan, 
Sadeghpour & Dalgarno values are based. Furthermore 
the fit formula is non-continuous at 85 eV. For energies 
between 30 and 85 eV we find that the cross section can 
be represented by 



CTbf,H2(a;) = y^^ajX 



(6) 



where x = E/15AeV and where crbf is given in Mbarn 
(= 10^^^ cm^). The fit coefficients Ui are given in Ta- 
bic Bl. The maximum deviation between the fit and the 
tabulated data is less than 1%. We note that the pho- 
toabsorption cross section of H2 is ~ 2.85(Th which results 
in an increase over previous estimates of ctism for regions 
where molecular hydrogen is an important contributor to 
the total absorptivity (i.e., below ~lkeV). Molecular ef- 
fects were not included in the earlier models, so essentially 
these models assumed cTbf.Ha = 2crbf,H- 

2.3. Abundances 

Naturally, the assumed elemental abundances are of cru- 
cial importance for the computation of (Tism • While many 
measurements of gas-phase abundances have been made, 
measurement of the total gas-plus-dust abundance of the 
ISM is very difficult. As a result, solar abundances, i.e., 
abundances determined from analysis of the solar photo- 
sphere or (carbonaceous) meteorites, have generally been 
used as the reference abundance for the ISM (Anders & 
Ebihara 1982; Grevesse & Anders 1989; Anders & Grevesse 
1989; ShuU 1993, and references therein). For reference, we 
list in Table B2 what we consider to be the best estimates 
of these "local" ISM abundances. We list the logarithmic 

^The abundances available in XSPEC are those of Anders & Grevesse (1989), Feldman (1992), and Anders & Ebihara (1982); in addition, 
user-specified abundances are possible. The BM92 model allows the user to vary the abundances with respect to these default abundances. 



abundances by number, normalized to a hydrogen abun- 
dance of log Ah = 12, based on the meteoritic abundances 
of Anders & Grevesse (1989). We updated the abundances 
of C and N using the results of Grevesse et al. (1991) and 
Grevesse & Noels (1993), respectively. The solar abun- 
dance of iron has been the subject of several discussions 
in recent years (Raassen & Uylings 1998; Kostik, Shchuk- 
ina & Rutten 1996; Biemont et al. 1991, and references 
therein). We adopted logApo = 7.50 ±0.05 as recently de- 
termined by Grevesse & Sauval (1999) in a critical rcx^valu- 
ation of the available data on Fe I lines. The uncertainty of 
the solar abundances in Table B2 is ~ 0.06 dex or smaller 
with the exception of log and logj4Ar, for which the 
uncertainty is 0.1 dex. 

As more abundance measurements have been made o\it- 
side our solar system, it has become apparent that the 
total gas plus dust ISM abundances are actually lower 
than the solar abundances (Sofia, CardcUi & Savage 1994; 
Savage & Scmbach 1996, and references therein). Stud- 
ies of the carbon abundance reviewed by Snow & Witt 
(1995) strongly indicate a subsolar abundance of carbon 
in the ISM, with a probable value of ~ 70% solar. Simi- 
lar reductions of 20% to 30% with respect to solar abun- 
dances are also indicated for the other metals (Snow & 
Witt 1996; Savage & Sembach 1996). This trend is consis- 
tent with new measurements of the gas-phase abundance 
of O, C, and N obtained with the Goddard High Reso- 
lution Spectrograph (GHRS) on the Hubble Space Tele- 
scope (Meyer, Jura &: Cardelli 1998; Meyer, Cardelli & 
Sofia 1997; CardeUi et al. 1996). 

Possible differences between solar and ISM abundances 
had been noted already at the time of publication of pre- 
vious models for ctism, and the previous computations of 
ciSM closely reflect the historical changes in the under- 
standing of these local abundances. Ride fc Walker (1977), 
for instance, used abundances that were greater than solar 
abundances in order to better represent abundance values 
in an ISM assumed to be chemically enriched by galac- 
tic evolution, and the model of BM92 was specifically de- 
signed to allow adjustments to the assumed abundances^. 
Now that a consistent picture of lower ISM abundance 
seems to have emerged, the cosmic abundance needs to be 
revised in models of ctism- To that end we adopt default 
recommended abundances for the ISM. These abundance 
values are listed in Table B2, and refer to the total abun- 
dance of the elements in either the gaseous or grain phase 
of the ISM. For C, N, and O, we use the values of Cardelli 
et al. (1996), Meyer, Cardelh & Sofia (1997), and Meyer, 
Jura & Cardelli (1998), respectively, while for all other ele- 
ments we use the abundances given by Snow & Witt (1996, 
Mg, Si, S, Ca, Ti, Cr, Fe, Ni), or, where no value was rec- 
ommended by these authors, we assume an abundance of 
70% with respect to our adopted solar abundances. 

We arc aware that these abimdances are still very much 
uncertain; the uncertainties in the general ISM abun- 
dances are still on the order of 0.1 dex or higher and also 
strongly depend on the line of sight. However, we feel that 
it has been adequately established that solar abundances 
arc enriched with respect to the ISM, and that for most 
astronomical work using the ISM abundances in Table B2 
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is preferable. Since one may anticipate the need for fur- 
ther revision of the interstellar abundances, we have also 
included an option to change these default abundances in 
the XSPEC implementation of our model (Appendix B), 
as BM92 did. 

Due to the appreciable difference in the X-ray cross sec- 
tions for molecular and neutral hydrogen (§2.2), it is nec- 
essary to include the contribution of molecular hydrogen. 
About half of the hydrogen in the Galaxy as a whole is 
molecular; however, it is not uniformly distributed (Shull 
& Beckwith 1982). Surveys of low galactic latitudes indi- 
cate that the ratio between H2 and H I is strongly depen- 
dent on the radial distance from the galactic center (Brinks 
1990; Bronfman et al. 1988). UV observations of H2 ab- 
sorption lines with Copernicus and with the ORFEUS II 
(orbiting and retrievable far- and extreme ultraviolet spec- 
trometer) echelle spectrograph indicate that about 20 . . . 
25% of the hydrogen is in H2 in the local ISM (Savage 
et al. 1977; Gringel et al. 2000). These results are very 
line-of-sight dependent, though, with some lines of sight 
being completely free of H2 . Because of these large fluctu- 
ations in the H2/H I ratio, it is thus difiicult to choose a 
"typical value" . In order to include the effect of molecular 
hydrogen on aisM we therefore choose to use the value of 
20% as the default value. To more accurately account for 
the opacity of molecular hydrogen, however, we strongly 
recommend that this ratio be set to that corresponding to 
the linc-of-sight to the source. The abundances of other 
molecules are much smaller so that their inclusion is only 
necessary for work along line of sights with very high Nu 
and hence not within the scope of this paper. 

2.4. Grain Physics 

Evidence for the presence of dust grains in the ISM 
comes from the observed interstellar extinction (Fitz- 
patrick 1999; Mathis 1990), as well as from the fact that 
the observed abundances of the elements in the ISM are 
generally less than the abundances assumed for the ISM 
(Savage & Sembach 1996). The ratio between the observed 
and the assumed abundances is called the "depletion", (3z, 
of the elements. It is generally assumed that the depletion 
indicates the presence of dust in the ISM, especially since 
the most depleted elements, C, O, and Si, are thought to 
be good candidates for the formation of solids in the ISM 
(Sofia, Cardelli & Savage 1994; Snow & Witt 1995; Mathis 
1996, and references therein). 

We have adopted a simplified grain model which as- 
sumes spherical grains in an MRN size distribution 
(Mathis, Rmnpl & Nordsieck 1977), 

dngr(a) _o r , ^ ^ 

— — — oc a where Omm < a < ctmax, (7) 

where typically grain size ranges from 0.025 /xm to 0.25 /xm 
(Draine & Lee 1984). While the grain composition and 
structure is still a matter of debate (Mathis 1996; Dwek 
1997; Smith & Dwek 1998; Wolff, Clayton & Gibson 1998), 
here, we assume composite ("fiuffy") grains, consisting of 
vacuum inclusions in solids. Thus, the grain density in our 
model (p = Igcm^^) is less than the p = 2gcm~'^ used 
in previous models (Fireman 1974; Ride & Walker 1977, , 
MM83). 

We assume chemically homogeneous grains consisting 
of silicates, graphites, and oxides. We list the adopted 



chemical composition of the grains in Table B2 in terms 
of the "depletion factor" 1 — the ratio between the gas 
abundance and the total interstellar abundance of the el- 
ements. Table B2 is based on the compilation of measure- 
ments given by Shull (1993), which is in general agreement 
with later GHRS measurements (Savage & Sembach 1996). 
Contrary to Shull, we do not assume N to be depleted into 
grains as is suggested by the observational evidence (Sav- 
age & Sembach 1996; Meyer, Cardelh & Sofia 1997). As 
recommended by Sofia, Cardelli & Savage (1994), we also 
assume a larger abundance of Fe in grains than that as- 
sumed by Shull (1993). For cobalt the values determined 
by Mullman et al. (1998) are used. 

To illustrate the large uncertainties associated with the 
grain composition, we also list in Table B2 the depletions 
used by Ride & Walker (1977) and MM83, the previous 
works which included grains in the computation of (tism- 
The depletion factors of Ride & Walker (1977), as shown 
in Table B2, are based on gas-abundance measurements 
using Copernicus data. Only one of the depletion factors 
of MM83 is based on observation, that for oxygen (de Boer 
1979, 1981); for the rest, MM83 considered elements to be 
either completely depleted into grains or entirely in gas 
form. This simplifying assumption represents the maxi- 
mum effect of grains on the ISM opacity. 

We note that grains are included here partly for com- 
pleteness and partly to indicate how changes to the current 
ISM grain model might affect observations made with the 
more recent X-ray observatories, and not as an attempt to 
model the solid state in the interstellar medium. While a 
"perfect" grain model meeting all abundance and obser- 
vational requirements remains elusive, our grain model at 
least adequately reproduces the observed extinction, emis- 
sion features, and depletions of the diffuse ISM (Mathis 
1996). For more detailed work, the assumption of spheri- 
cal grains must be modified (Mathis 1990). 

3. CONCLUSIONS AND SUMMARY 
3.1. Results 

Using equation (1) and the assumptions described 
above, we compute cigM. The results are shown in Fig. 1, 
where we display the X-ray absorptivity as ctism^'^ to em- 
phasize deviations from the proportionality of a^f- 
By reading the plot as a "bar-diagram" , it is easy to esti- 
mate the relative importance of the contributors to ctism- 
For energies above the oxygen K-edge at ~0.5keV(see 
Tab. B2), the X-ray opacity is dominated by the metals, 
and H and He are relatively unimportant. Below IkeV, 
C, N, O, and Ne are the important absorbers, while above 
1 keV, Si, S, and Fe are important. 

Also clear from Fig. 1 is that the effect of grains on 
Ctism is small for a standard MRN distribution. The ef- 
fect of grains is found to be less than previous estimates 
such as MM83 in part because our grain model consists of 
lower-density porous grains and in part because we calcu- 
late CTgrains for an MRN distribution of grain sizes rather 
than by choosing an average grain size such as 0.3 /im 
as used by Fireman and others. Thus the greatest num- 
ber of grains lies in very small grains with r ^ 1 above 
~ IkeV, so that self-shielding is not important in these 
grains. Consequently, for an MRN distribution there are 
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0.1 1.0 10.0 

E [keV] 

Fig. 1. — Absorptivity per hydrogen atom of the ISM using the assumptions described in the text. The dotted line is the absorptivity 
including grains with a MRN distribution, and the dashed line is the absorptivity assuming that all grains are of radius a = 0.3 /im. The 
inset shows the cross section without the multiplication with E^. We also illustrate the contribution of hydrogen and hydrogen plus helium 
to the total cross section. The contribution of the H2 cross section to the total hydrogen cross section is indicated by the dot-dashed line. 




Fig. 2. — Absorptivity per hydrogen atom of the ISM for our adopted ISM abundances (solid line) and for the solar abundances of Tab. B2 
(dashed line). For clarity, the cross section has been multiplied with E^; however, the 3/-axis is logarithmic here, not linear as in Fig. 1. The 
revised abundances produce changes in ctism of up to 30% with respect to the solar abundances. 
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Fig. 3. — Optical depth of grains for a hydrogen column ATjj = 1.0 x 10^" cm~^ (see Appendix A), and for grains of radius a = 0.25 jxra 
(dot dashed line) and a = 0.025 ixm. (dashed line) as a function of X-ray energy. The total mass in grains is constant. The shielding reduces 
the optical depth at low energies and for large grain sizes. The labels denote the K-edges of the elements indicated. 



only small differences between the optical depth of an en- 
tirely gas-phase ISM and the optical depth for an ISM 
in which some of the gas atoms have been depleted into 
grains. This is illustrated in Fig. 3, where we plot the grain 
optical depth as a function of energy for grains of radius 
0.25 /v,m and 0.025 /xm, as well as for an MRN distribution 
of grains. These optical depths are calculated using Ap- 
pendix A and considering a hydrogen column density A^h 
of 1.0 X 10^'^cm^^. The total mass in grains is the same for 
all three grain models shown. Clearly, the self-blanketing 
factor affects the opacity more for low energies and larger 
grain sizes. Since the change in grain optical depth for 
the MRN distribution is quite small, grains only slightly 
modify the absorptivity of the ISM. 

We stress, however, that the grain model as used here 
applies to the diffuse ISM rather than to dense regions such 
as inner molecular clouds or young stellar objects, where 
grains might play a much more important role (Mathis 
1990). By changing the range of grain sizes from the range 
adopted here for the diffuse ISM to a grain distribution 
with larger grains, our formalism (including the XSPEC 
routine described in Appendix B) can be used to treat 
these regions as well. As an example, we show in Fig. 1 
the influence of grains with a = 0.3 /im. In this case. X-ray 
observations will allow the detection of the grains due to 
the large influence of self-shielding, even with instruments 
that do not have the high energy resolution to see the 
solid-state resonance effects close to the absorption edges. 

3.2. Comparison with previous models 



As we outlined in Sect. 2.2 and 2.3, the major differences 
between this model and the previous models are the up- 
dates to the atomic cross sections and the assumed abun- 
dances. With the use of simulated observations, we can 
examine the impact of an updated (Tism on the interpre- 
tation of observations. 

In order to demonstrate the effect of the updates to the 
cross sections, we compare our values of (Tism with those 
of MM83 and BM92. In Fig. 4a we show the relative de- 
viation between our formulation and that of the earlier 
models, using our adopted cross-sections and the assumed 
abundances of the MM83 and BM92 models, respectively. 
Thus, the figure shows the cumulative effect of the up- 
dated cross sections. The differences betweciii the models 
at photon energies below 0.3 keV are mainly a result of 
the uncertainty in the He cross section. Additional peaks 
at the K- and L-shell edges are due to ^ lOeV discrep- 
ancies in the assumed edge energy between these models 
(§2.2) and are not relevant for practical work with moder- 
ate resolution instruments. The difference in the values of 
fTbf produces the curvature of the relative deviation. The 
missing subshell absorption edges in the MM83 model are 
also readily apparent. 

In Fig. 4b, we again show the relative deviation between 
our model and that of the earlier models, but this time us- 
ing our default solar abundance (Table B2) for the abun- 
dances in our model, including the additional elements not 
contained in the MM83 and BM92 models, as well as the 
II2 molecule (important especially below ^ 300 eV). In this 
figure, all models are with "solar" abundances, but the dis- 
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Fig. 4. — Relative deviation of the cross sections of BM92 (solid line) and MM83 (dashed) from the model described here. The relative 
computation is given as {ai — a)/cr where i =BM92 or MM83, and where cr is our formulation. The spikes at the K- and L-absorption 
edges are due to slight differences in the assumed edge energies, a: Relative deviation computed using the abundances of BM92 and MM83, 
respectively. This figure highlights the effect of the change in the cross-sections, b: Relative deviation to our model, computed using the solar 
abundances of Tab. B2, illustrating the cumulative effect of using different default solar abundances and changing cross sections. See §2.2 for 
further discussion. 



8 



Wilms, Allen, & McCray 



agreement between the "solar" abundances used by the dif- 
ferent authors results in 10% of the difference between the 
three models compared in Fig. 4. However, the difference 
resulting from different values of the "solar" abundance is 
small compared to the difference that results from using 
our default abundance values which are reduced with re- 
spect to solar abundances. The discrepancy in ctism which 
results from the use of our default ISM abundances instead 
of the solar abundances indicated in Table B2 is shown in 
Fig. 2. 

In order to illustrate the observable changes introduced 
by our updated model of ctism, we studied these effects us- 
ing simulations of observations. We simulated the obser- 
vation of a bright X-ray source which has a pure power-law 
spectrum with photon index T = —1.7, and a normaliza- 
tion of ApL = 2 phcm~^ s~^ keV~^ at 1 keV. Such a source 
spectrum is characteristic, e.g., of Cygnus X-1 in the hard 
state (Dove et al. 1998). Wc then used two values of hydro- 
gen column, A^H = 6 x 10^^ cm~^ and Nu = 6 x 10^^ cm~^, 
to produce the absorbed spectrum using our model for 
ciSM, using our assumed interstellar abundances. The two 
values for A^h are necessary to illustrate the effects of A^h 
both in the low energy band (small Ajj) and in the high 
energy band (large A^h)- We simulated an observation by 
folding the photon spectrum through the detector response 
matrix and adding Poisson noise. To ensure comparable 
signal-to-noise ratios in all simulated observations, the ex- 
posure times were chosen such that 10^ photons were con- 
tained in the resulting spectrum. The detectors we chose 
to represent current instruments were SIS on the Advanced 
Satellite for Cosmology and Astronomy (ASCA) for the 
low A^H simulation, and the Proportional Counter Array 
(PCA) on the Rossi X-ray Timing Explorer (RXTE) for 
the large A'h simulation. The detector chosen to represent 
the quality and energy resolution of the future missions is 
the European Photon Imaging Camera (EPIC) pn CCD on 
XMM-Newton. In neither case did we take photon pileup 
or other instrumental effects such as a saturation of the 
telemetry bandwidth into account. 

The resulting simulated photon spectrum was rebinned 
and analyzed using the MM83 and the BM92 models for 
(TisM- In Fig. 5 we present the ratio between the data 
and the resulting best-fit model using the MM83 model 
(the BM92 model without varying the abundances gave 
qualitatively the same results). For the PCA and the SIS, 
the fits with the MM83 model gave an acceptable descrip- 
tion of the data. Our low A^h simulation shows that the 
ASCA SIS is able to detect the discrepancy in the oxygen 
abundance between the solar and the ISM values. The 
other abundance differences are not measurable in reason- 
able amounts of observing time. We note in passing that 
Ebisawa et al. (1996) indeed found residuals in an observa- 
tion of the black hole candidate Cygnus X-1 that are very 
similar to those shown in Fig. 5, which points towards 
an oxygen abundance of ~80% solar in the direction to 
Cyg X-1. 

For the large A^h simulation, the fit to the simulated 
RXTE PCA spectrum recovered the input parameters for 
r and for ApL to two significant digits and found a hydro- 
gen column of A^h = 4 x 10^^ cm~^. The residuals shown 
in Fig. 5 show some structure below SkeV that is caused 
by underestimating A^h, which results in the non-optimal 



value of the fit (x^ed = ^ ■''')■ However, the PCA has 
significant systematic uncertainties in this region (Wilms 
et al. 1999), and the deviation seen in our simulation would 
probably go undetected in a real observational situation. 

Although we have not done an exhaustive search of the 
whole A^H space, the above result convinces us that for 
most applications using the current moderate resolution 
instruments, the models of MM83 and BM92 are suffi- 
cient. The physical cause for the underestimate of A^h is 
that for energies above 1 kcV, ctism is dominated by the 
metals. Since the abundance ratios of the metals are a 
result of stellar nucleosynthesis, reducing the metal abun- 
dance with respect to hydrogen from the MM83 or BM92 
values to our values can be compensated by assuming a 
different (lower) hydrogen column. Nor does the slightly 
different shape of abf that results from the updated cross 
sections produce a noticeable effect in the PCA simula- 
tions; the energy resolution of the PCA is insensitive to a 
difference of this magnitude. 

Our simulations for the EPIC pn camera, on the other 
hand, show deviations on the 5% level when fitted with 
the MM83 and BM92 models (Fig. 5). These changes are 
clearly due to the update of the abundances and cross sec- 
tions. It is for instruments such as the EPIC pn camera 
that we recommend that our model be used. The EPIC 
observation simulated here is quite short (the raw EPIC 
pn count-rate of Cygnus X-1 will be ~ 2400 cps, so that 
the spectrum used for modeling here has an exposure time 
of only 4 ksec - shorter than the minimum XMM-Newton 
observing time of 5 ksec), and the ratio plots in Fig. 5 al- 
ready show the characteristic curvature due to the devia- 
tion in the shape of the cross section assumed in the MM83 
models. We note that a similar behavior is also present 
when fitting the simulated spectrum with the BM92 data. 
When using the BM92 cross sections, however, it is possi- 
ble to obtain a good description of the overall absorption 
by treating the abundances as free parameters. In this 
case the "wavy" residual structure still remains although 
the relative abundances found for the elements do come 
close to the abundances that we assumed in our modeling 
of (JisM- As a result, wc recommend using our model when- 
ever high precision work on the l%-2% level is necessary, 
as well as in cases when the absorbing column resulting 
from the spectral fitting process is to be physically inter- 
preted (i.e., when comparing X-ray A^h values with radio 
data, or when modeling an observed X-ray spectrum for 
which the A^h value is known by other means). 

3.3. Caveats 

While wc feel that the model for (Tism as presented in 
the previous sections is well-suited for most X-ray astro- 
nomical work, several astronomical or physical issues had 
to be neglected in the modeling process. In particular, our 
model does not include the effects of the warm phase or of 
the ionized phase of the ISM. In addition, wc ignore Thom- 
son scattering, which becomes important at moderate to 
high Thomson optical depths and E >4 keV. These effects 
have not been included largely because they are negligible 
for most situations in which our model is likely to be used. 
Furthermore, the model assumes that there is no varia- 
tion of the abundances along the line of sight. Finally, we 
note that energy-dependent scattering by dust grains has 
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Fig. 5. — Ratio of simulated X-ray data absorbed with our model for ajsM to the best fit model using a power-law model and the MM83 
models, for the ASCA SIS (triangles), RXTE PCA (diamonds), and the XMM-Newton EPIC pn CCDs. The ASCA SIS simulation with our 
model is shown for energies below 2 keV. In all cases, the error bars represent Poisson noise. While older models produce adequate fits for the 
current instruments (with the exception at the O edge, see text), differences of ~ 5 . . . 10% are found with the EPIC pn simulations. 
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not been considered here since an XSPEC routine (dust) 
already exists to treat it. 

The warm and ionized phases of the ISM have been 
neglected here for a number of reasons. Large parts of 
the ISM arc assumed to be (moderately) ionized, and a 
good model for aisM should in principle include the ion- 
ized phase. However, there is even less agreement on the 
ionization state of the ISM as a whole than there is on 
the grain phase, which was discussed in §2.4. The major 
reason for this disagreement is a good one: the ioniza- 
tion of the ISM depends strongly on the line of sight, so 
that any general model for the ionization structure of the 
galaxy is prone to fail. It appears, however, that most of 
the warm phase of the ISM is only moderately ionized, so 
that our model will still be applicable except for photon 
energies close to the K-edges. For oxygen, e.g., the K-edge 
energy increases from 0.54 keV for O I to 0.61 keVfor O IV, 
and the photoionization cross section is also dependent on 
the ionization stage. It should be noted, though, that the 
ionization of hydrogen will strongly influence cism for en- 
ergies below O.lkeV. For instance, in their study of soft 
X-ray emission from other galaxies, (Cui et al. 1996), as- 
sumed the ratio between H I and H II to be 1. Using the 
code presented in Appendix B, this effect may be taken 
into account by changing the relative abundances of He 
and the metals with respect to hydrogen. 

Similar line of sight effects hold also true for the abun- 
dances. In the galaxy, measurement of abundanc;cs in HII 
regions and planetary nebulae has suggested the existence 
of a radial abundance gradient, particularly at small galac- 
tic radii (Rana 1991; Kaufcr et al. 1994, and references 
therein). When computing the absorption of any back- 
ground object through such a medimn with varying abun- 
dances, the abundances measured in the X-rays will be the 
average of the abundance variations along the line of sight. 
Since the abundances in our model are allowed to vary, it 
is possible to take this effect into account. Obviously, in 
such a case the assumed abundance set with respect to 
which the relative abundances are measured is a matter 
of taste, and using the solar abundances of Tab. B2 as a 
baseline instead of the ISM abundances might be a better 
choice. 

Our model also ignores the influence of Thomson scat- 
tering because for many observations it will not make a 
substantial contribution. This is not to say that it is al- 
ways negligible; Thomson scattering is the dominant phys- 
ical process for the attenuation of X-rays above ~10keV. 
As a rough guideline, Thomson scattering should always 
be included in the modeling process when TVh > 10^^ cm~^ 
(see, e.g., Stelzer et al. 1999). Detailed instructions based 
on Monte Carlo modeling of Thomson scattering in mod- 
erately optically thick media are given by Yaqoob (1997). 
Note, however, that Yaqoob assumes the electrons to be 
free. When the electrons are boimd, it is more appropriate 
to consider incoherent scattering, i.e., Compton scattering, 
off the bound electrons. The cross sections for incoherent 
scattering and the Klein-Nishina cross section are slightly 
different (on the percent level), which should be taken into 
account. Analytical estimates for this effect have been 
given by Gorshkov, Mikhailov & Sherman (1973) for the 
H atom; for other elements tables are available (Hubbell 
et al. 1975; Henke, GuUikson & Davis 1993). For most 



work, however, these corrections are negligible. 

Finally, we also stress that our model, like its prede- 
cessors MM83 and BM92, is not well suited for extremely 
high-precision work, e.g., high signal-to-noise spectra with 
the grating spectrometers on XMM- Newton and Chandra. 
The energy resolution of these very-high-resolution instru- 
ments is so high that the X-ray absorption fine structure 
(XAFS) within grains and molecules might become rele- 
vant, as well as resonance effects close to the absorption 
edges in the absorption cross section itself. For the inclu- 
sion of the resonances much higher resolution cross sec- 
tions than those available here are required. XAFS effects 
become noticeable especially in regions of large TVh where 
grains are an important ingredient. For these studies, ded- 
icated modeling codes are being devised, which should be 
used instead. See Forrey, Woo & Cho (1998), Woo, Forrey 
& Cho (1997), Woo (1995), and the references therein, for 
further details. 

3.4. Summary and Outlook 

We have presented an improved model for the X-ray ab- 
sorption in the interstellar medium valid above ~ 100 eV 
by updating the relevant photoionization cross sections 
and by including the effect of the H2 molecule and a re- 
vised set of recommended elemental abundances for the 
interstellar medium. Since the abundances of the metals 
are smaller than previously estimated, our values of ctism 
are generally smaller than those used previously. We have 
shown that these updates will be relevant for instruments 
to be available as this paper goes into print. We have also 
reconsidered the effect of grains in the ISM by improving 
the grain model and the calculation of Tgrains) taking into 
account recent studies of the physical properties of inter- 
stellar dust grains. Our results show that the change in 
opacity due to grains on cjism for realistic grain models and 
realistic chemical composition of the grains is small. As 
described in Appendix B, our calculation of (Tism is avail- 
able as a subroutine that can be used with X-ray data 
analysis packages such as XSPEC. 

This work highlights the point made by Shull (1993) 
and others, that great uncertainties remain in much of the 
atomic physics data needed to interpret X-ray astronom- 
ical observations. The effective areas and energy resolu- 
tions of instruments recently launched will be more sensi- 
tive than the available atomic data for the X-ray energy 
range can accommodate, and the situation is even worse 
for possible future instruments such as Constellation-X 
or XEUS. The differences between the currently available 
cross sections discussed in §2.2 are still at the 5% level. We 
hope that the advent of the new instruments will also initi- 
ate new attempts to more accurately determine the atomic 
data on wliic;li X-ray astronomy can rely. Wo will there- 
fore continue over the next years to update the database of 
the computations presented here and periodically release 
new versions of the code presented in Appendix B that 
will incorporate the improvements in the abundances or 
the atomic physics. 
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APPENDIX 



SELF-SHIELDING IN THE GRAINS 



To derive the optical depth of the grains (eq. [5]) we consider a population of grains of homogeneous chemical composition 
with size a, geometrical cross section Cgeom, column density (N), and average photoionization cross section (a) (per atom). 
For grains of realistic shape, {N) is a complicated function of a. Since the influence of grains on ctism is small, however, 
for simplicity we can assume that the grains are spheres with radius a so that 
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where p is the mass density of the grain, n the number density of the grain, and ji the mean molecular weight of the 
grains constituents. Using the notation of Section 2.1, fj, is given by 

lJi = Y,Az-Pz-l^z /Y.^z-I3z (A2) 
z ' z 

Here, fiz is the molecular weight of element Z (Tab. B2). For the default abundances and depletions of Tab. B2, 
fi = 18.2 amu = 3.0 x 10~^^g/atom. Assuming p = Igcm^-^, the grain number density is n = 3.3 x 10^^ atoms cm~^. 
Ignoring solid state effects, the photoionization cross section of the grain material is given by 

(a) =J2^z-0z- <Ju{Z) lY^Az-fiz (A3) 

z ' z 

For a column density iVgrains of grains of a single size along the line of sight, the optical depth for photoabsorption of the 
grains is the product between the geometrical optical depth of the grains, A^grainsfgeom, and the fraction of photons that 
encounter the grain and are absorbed: 

rains = A^grainsO-goom (1 " exp (-(cr)(7V))) . (A4) 



To 



In computing ctism for grains of a single size it is sometimes convenient to express Tgrains in terms of the optical depth 
that would be observed if the constituents of the grain were in the gaseous phase, Tgrains = /-^gas(o'). Hence, 

-^gas(cr) W){N) 

since the column density of the grain is Afgas/(AfgrainsO'geom)- This factor is called the self-blanketing factor (Fireman 
1974). 

For the general case that the grains have a size distribution dngr(a)/do, where 

'^'''''da=l (A6) 



/ 

Jo 



(A9) 



da 

a similar derivation shows that 

Tgrains - A^nCg ^^^Tgeom (1 " cxp (- (a) (iV) ) ) da, (A7) 
where ^g is the number of grains per hydrogen atom along the line of sight. For spherical grains, 

ig = Y.^z-(iz-p^zl J^-p-^na'-da (A8) 

For a power-law size distribution, dngr(a)/da = ka~P, where amin < a < aamax, we can express Tgrains in terms of special 
functions. Inserting Cgeom = Tra^ into equation (A7) and performing a partial integration, we obtain for p < 4 

r ■ - NT,f / ^ ( 1 - p-4(<T>na„i„/3'\ 3_ (l- e-4(<7>na„.ax/3'\ 

Tgrams — JVH^g „i p.gl-L e I 3 11 6 1 

+ {A{a)n/ 3)f-3 (4 _ 4(a)na„ax/3) - 7 (4 - 4(a)na„in/3)] }, 
where the incomplete gamma function is given by (Abramowitz & Stegun 1964, section 6.5) 

7(a,a:)= / e-*t"-Mi (AlO) 
Jo 

and can be evaluated using algorithms presented by Press et al. (1992, Section 6.2). 

IMPLEMENTATION INTO XSPEC 

Using the formalism and data described in Section 2, wc have written subroutines that can be used in conjunction with 
the popular X-ray data analysis package XSPEC (Arnaud 1996). These subroutines are available at: 
http : //astro . uni-tuebingen . de/nh/ 

Given the uncertainty of many of the parameters entering the model, there are several interfaces to the subroutine. 
The main spectral model, tbabs, is intended as a replacement for the earlier wabs and phabs models, as well as their 
redshifted versions. This model assumes the abundances and depletion factors of Table B2. A subroutine ztbabs for 
modeling redshifted absorption is also available. The latter subroutine does not account for the dust component, as X-ray 
observations with a measurable dust influence will most probably only include observations of galactic objects, at least 
within the next decade. 

In addition, the XSPEC model tbvarabs provides a full interface to the X-ray absorption model presented in Section 2 
and allows the user to change all abundances, depletion factors, and grain properties. We stress that this model is very 
powerful and therefore prone to misuse. Caution is required when fitting data with too many free parameters. 



X-Ray Absorption in the ISM 



Table B1 

Fit coefficients Ui to equation (6) for the H2 photoabsorption cross section. 



ao 


ai 


a2 


as 


04 


05 


0.664 


-11.768 


78.118 


-231.339 


368.053 


-189.953 



Table B2 

Molecular weight, K edge energies, abundances and depletion factors 1 — f3z for the abundant elements. 



Element 






I2 + I0 


gAz 












keV 


solar'^ 


ISM'i 


this paper 


MM83*' 


Ride77s 


1 H 


1 




12.00 


12.00 


1.0 


1.0 


1.0 


2 He 


4 




10.99 


10.99 


1.0 


1.0 


1.0 


6 C 


12 


0.29 


8.60 


8.38 


0.5 


0.0 


0.2 


7 N 


14 


0.41 


7.97 


7.88 


1.0 


0.0 


0.5 


8 


16 


0.54 


8.93 


8.69 


0.6 


0.75 


0.5 


10 Ne 


20 


0.87 


8.09 


7.94 


1.0 


1.0 


1.0 


11 Na 


23 


1.08 


6.31 


6.16 


0.25 


0.0 




12 Mg 


24 


1.31 


7.59 


7.40 


0.2 


0.0 


0.2 


13 Al 


27 


1.57 


6.48 


6.33 


0.02 


0.0 




14 Si 


28 


1.85 


7.55 


7.27 


0.1 


0.0 


0.5 


15 P 


31 


2.15 


5.57 


5.42 


0.6 






16 S 


32 


2.48 


7.27 


7.09 


0.6 


0.0 


0.7 


17 CI 


35 


2.83 


5.27 


5.12 


0.5 


0.0 




18 Ar 


40 


3.20 


6.56 


6.41 


1.0 


1.0 


0.5 


20 Ca 


40 


4.04 


6.34 


6.20 


0.003 


0.0 




22 Ti 


48 


4.97 


4.93 


4.81 


0.002 






24 Cr 


52 


5.97 


5.68 


5.51 


0.03 


0.0 




25 Mn 


55 


6.55 


5.53 


5.34 


0.07 






26 Fe 


56 


7.12 


7.50 


7.43 


0.3 


0.0 


0.2 


27 Co 


59 


7.73 


4.92 


4.92 


0.05 






28 Ni 


59 


8.35 


6.25 


6.05 


0.04 


0.0 





^molecular weight in amu, 

edge energy (rounded after Verner & Yakovlev 1995) 

'^solar abundance (Anders & Grevesse 1989; Grevesse et al. 1991; 
Grevesse & Noels 1993, see also §2.3), 

''adopted abundance of the ISM based on Snow & Witt (1996), Cardelli 
et al. (1996), Meyer, Cardelh & Sofia (1997) and Meyer, Jura & Cardelh 
(1998), see §2.3; 

**ratio of gas abundance to total ISM abundance, 1 — Pz, using our 
adopted abundances (§2.4). For comparison, the older values of 

f Morrison & McCammon (1983), and 

sRide & Walker (1977) are also hsted. 



